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The origin of the spin gap in the underdoped cuprate superconductors is still mys- 
terious. Experimental evidence from neutron scattering and NMR experiments in- 
dicates that the spin gap might be present only in the bilayer compounds. A naive 
calculation for a two plane Heisenberg model locates the order-disorder transition 
only for very large exchange coupling between the bilayers. We propose a interplane 
relaxation experiment which might quantitatively estimate the strength of the bilayer 
coupling. We make detailed predictions for the size and the temperature dependence 
of the interplane relaxation rate. 



1. INTRODUCTION 

Recently the low energy spin dynamics of the high T c superconductors has 
attracted a lot of interest. Some of the theories for the high temperature super- 
conductors are based on low frequency spin fluctuations. 1 The unusually strong 
temperature dependence of the measured spin susceptibility, Xs(T) of some of these 
compounds (YBa2Cu3C>6 63, YBa2Cu408, • • •) rules out a Fermi liquid like ground 
state with a temperature independent spin susceptibility. The decreasing spin sus- 
ceptibility indicates the closeness of a critical point. Many people have proposed 
that this critical point should be identified with the order-disorder transition in the 
two dimensional Heisenberg Quantum Antiferromagnet. 2 The reason why these sys- 
tems are close to this transition is however still under discussion. Millis and Monien 3 
proposed that the bilayer coupling is driving the system towards the quantum dis- 
ordered phase whereas Sokol and Pines 4 proposed that doping alone is sufficient to 
reach the quantum disordered phase. This critical point has recently been studied 
in great detail by Sachdev et al. and others. 5-8 Millis and Monien concluded from 
the analysis of existing experimental NMR data that the spin gap only appears in 
systems containing CuO'i bilayers. An additional hint that bilayer correlations are 
important in YBCO comes from the neutron scattering experiments by Tranquada 
et al. 9 which show nearly complete antiferromagnetic correlations of the spins in 
the bilayers. All this points in the direction of a sizable exchange coupling between 
the planes. The study of a model of two antiferromagnetically coupled Heisenberg 



planes reveals that the exchange coupling between the planes, Jj_ has to be of the 
order 2.5 x J, 7,8,10 the antiferromagnetic exchange coupling in the plane, which is of 
the order of J ~ 1000A'. This exchange coupling seems too large for a realistic ma- 
terial. Several groups 11-13 have attempted to clarify mechanisms which could lead 
to an enhancement of the singlet correlations between the bilayers. In this paper we 
analyze theoretically in more detail a recently proposed experiment, 14 which makes 
it possible to directly measure the transition of the "correlated plane regime" to the 
"uncorrelated plane regime" and to obtain an estimate for the exchange coupling 
between the planes Jj_ . 

2. INTERPLANE RELAXATION 

The basic feature of all CuO superconductors are the CuO planes. In 
YBa2Cu306+5 two layers form a so called bilayer which is separated from the next 
bilayer by a large distance. Each of the layers of the bilayers is completely equiva- 
lent in YBa2Cu306+5- Recently a new CuO superconductor has been synthesized 
in which the CuO'i layers of the bilayers are not equivalent. The unit cell con- 
tains two building blocks one of which has one chain, similar to YBa2Cus07 and 
the other two chains, like YBa2Cu408- The structure is sketched in Fig. 1. The 




Fig. 1. Schematic structure of Y2Ba4Cu70i5 

important point is that each of the CuO'i layers are attached to a different chemi- 
cal environment but are otherwise left intact. One plane of the bilayer is attached 
to the single chain block the other one to the double chain block. Some of the 
properties of this material have been studied by Stern et al. 15 In particular they 
find that the nuclear quadrupole resonance (NQR) frequencies of the Cu nuclear 
spins in the individual planes of the bilayers are different from each other. This 
allows one to access the spins in both Cu02 planes seperately. It is known that the 
transverse relaxation rate, I/T2, is dominated by an indirect exchange coupling via 
the conduction electrons. 16 A Cu nuclear spin at site i in plane (a), couples to an 



electron spin in the same plane, (a), via the hyperfine Hamiltonian 
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where A aa , a G x,y,z, is the on-site hyperfine coupling for the Cu 2+ spin and B 
is the isotropic transferred hyperfine coupling. The sum over j extends over the 
four neighboring Cu sites in the same plane. The nuclear Cu spin polarizes the 
neighboring electronic spins in the same plane which in turn produce an additional 
field for the neighboring Cu nuclear spins. This process generates a nuclear-spin 
nuclear-spin interaction of the form 
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V aa (Ri — Rj) is defined by its Fourier transform: 

V aa (q) = F aa ( q )Ka(Q)xi a ' a) (q) 



(3) 



where F aa (q) is the Fourier transform of the hyperfine coupling constant in Eq. 
(1) and X^ a ' a \l) is the static spin susceptibility for a single plane. For a Cu 2+ 
spin in a d x 2_ y 2 orbital the hyperfine interaction is anisotropic and largest in the 
c direction perpendicular to the planes. The Cu nuclear spin interaction, Eq. (2), 
dominates the dipol-dipol interaction of the nuclear spins because the spin suscepti- 
bility, Xs(l), is strongly enhanced at the antiferromagnetic wavevector in the plane. 
The interaction is giving the main contribution to the dephasing time of the Cu 
nuclear spins in one CuO'j layer. Pennington and Slichter 16 derived an expression 
for the dephasing time T2 for the case that the hyperfine coupling in the c direction 
is much larger than any other direction: 
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where n m is the density of the active NMR nuclei. The second term cancels the on- 
site contribution of the first term. This time has been measured in YBa2Cus07 and 
is in reasonable quantitative agreement with simple phenomenological theories. 16, 17 
In the special bilayer system discussed above it is possible to excite a spin in one 
plane and measure the response in the other plane. Then one can use Eq. (4) to 

define a interplane-relaxation time for the bilayer system. In Eq. (4) Xq "'^ nas to be 
(12) 

replaced by Xq where the spins are sitting in different planes. In this case we do 
not have to subtract the on-site term. To discuss the the interplane relaxation time 
we use the two plane Heisenberg model which exhibits a order-disorder transition 
as a function of the between plane coupling Jj_ . 

Many physical properties of the two plane Heisenberg model in the disordered 
phase can be calculated in an approximate way with the Schwinger-boson mean-field 
theory. Schwinger-boson mean-field (SBMFT) does not reproduce the numerical 
value of the critical exchange coupling between the planes, Jj_ (numerical value 



Jl ps 2.5J 8 ' 10 ' 18 vs. J ± ps 4.5 J in the SBMFT 19 ) but it does correctly give the 
main features of the phase diagram. The spin-correlation functions and a number 
of physical properties of the two plane Heisenberg model have been obtained in this 
approach. Here we are interested in the calculating the correlation function of a spin 
in one plane, with a spin in the second plane, S (2 'K The SBMFT is rotationally 
invariant in spin space. We therefore can confine ourselves to calculating the (S z S z ) 
correlation function. An external field in the z direction S z is coupling to the spin, 
S\ a \ in plane (a), site i via the Hamiltonian 

AH = J2hM a) (5) 

ia 

The elementary excitations of the system are the optical and the acoustic spin- 
waves. It is convenient to rewrite the Eq.(5) in terms of these excitations: 

AH = Y<{1 ~ <>\ + \ + <>\ } (6) 

where O q and O q are operators creating spin fluctuations symmetric and antisym- 
metric under interchange of the planes respectively. The only nonvanishing correla- 
tion functions are < O q O s _ q > and < O a q O a _ q >. In linear response theory we obtain 
for the spin-correlation function when both spins are sitting in the same plane: 

(Sf) S*}?) (q,u) = \ (x ss (c l ^)+X aa (^)) (7) 
where x aa { f li L0 ) an d X 88 ^!^) are defined as 

{• CO 
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and the index (1) is referring to the plane one. The spin-correlation function for 
both spins sitting in plane (2) is of course identical. 

For the spin-correlation function with the two spins in different planes we ob- 
tain: 

(S*U S^)( q ,u 1 ) = ^( X ss ( q ,u 1 )- X a % q ,u 1 )) (10) 

Now we are in the position to compare the interplane relaxation time, T21 
with the transverse spin relaxation time T2. At low temperatures the planes are 
coupled and only the acoustic mode is populated. In this regime x" is much larger 
than x aa ■ Comparing Eqs. (7) and (10) and we see that T2 and T^l are identical. 
With increasing temperature the optical mode is more and more populated until, 
at approximately T ps Jj_ the planes are decoupling and the optical and acoustic 
mode merge. In this temperature regime x aa an d X'" are equal. 

This transition is of course an artifact of the SBMFT and will be a interpla- 
neover if fluctuations are taken into account. For this reason we have studied the 
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Fig. 2. Temperature dependence of the interplane relaxation rate I/T21 for Jj_/J = 
2.5 on a 2 x 32 x 32 lattice and the ratio of of the interplane relaxation time, T2j_,to 
the transverse nuclear relaxation time, T2. 



two plane Heisenberg model using a recent improvement of the Hanscomb Quan- 
tum Monte Carlo method which is particularly suited for non frustrated systems. 20 
We have calculated the interplane relaxation time for various values of Jj_. Fig. 2 
shows the interplane relaxation time as a function of temperature for a Jj_ = 2. 5 J, 
close to the order-disorder transition. In this calculation we have used the values 
for the hyperfine coupling as determined from other NMR experiments. The in- 
terplane relaxation rate is strongly decreasing with increasing temperature. There 
is no strong feature visible at the temperature of T ps 2.5 J where the Schwinger- 
boson mean-field theory would predict the transition from the coupled plane to the 
decoupled plane regime. This has several reasons. The acoustic mode is becoming 
broad with increasing temperature and starts to overlap with the optic mode. At 
a temperature of 2.5 J it is very natural not to see a sharp transition. One might 
hope that for smaller values of Jj_ the crossover is more pronounced. The value 
of the interplane relaxation rate, I/T21, for a Jj_/J ~ 2.5 is much larger than the 
preliminary experimental findings by Mali and Stern. 21 We are currently exploring 
smaller values of Jj_ . The correlations between the planes seem to be strong even 
for rather small values of Jj_ . 



3. CONCLUSIONS 

We have analyzed an interesting experiment which might be able to pin down 
the value of the bilayer coupling in YBCO. The preliminary results of Mali and 
Stern 21 seem to indicate that T'j±_ is of the order of T%. Our recent calculation show 
that only a modest value of Jj_ is required to obtain a large T21. The values for Jj_ 
are in agreement with earlier estimates by 0. K. Andersen 22 from a bandstructure 
analysis. An experimental and theoretical comparison between T2 and T^l might 
help to identify the decoupling transition. 
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